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Abstract 

Measurements of the time-dependent CP asymmetry in the decay B^(t) — > 7r + 7r~ and its charge 
conjugate by the BELLE and BABAR collaborations currently yield = —0.46 ± 0.13 and = 

—0.74 ± 0.16, characterizing the direct and mixing-induced CP asymmetries, respectively. We study 
the implication of these measurements on the CKM phenomenology taking into account the available 
information in the quark mixing sector. Our analysis leads to the results that the ratio \P C /T C \ involving 
the QCD-penguin and tree amplitudes and the related strong phase difference 8 C = 5^ — 5^ in the 
B^/B^ — > 7r + 7r~ decays are quite substantial. Using the isospin symmetry to constrain |P C /T C | and 
cos(2(9), where 29 parameterizes the penguin- induced contribution, we present a fit of the current data 
including the measurements of S+~ and C+~. Our best-fits yield: a = 92°, (3 = 24°, 7 = 64°, \P C /T C \ = 
0.77, and 8 C = -43°. At 68% C.L., the ranges are: 81° < a < 103°, 21.9° < j3 < 25.5°, 54° < 7 < 75°, 
0.43 < \P C /T C \ < 1.35 and —64° < S c < —29°. Currently en vogue dynamical approaches to estimate the 
hadronic matrix elements in B — » txtx decays do not provide a good fit of the current data. 
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1 Introduction 



Precise measurement of CP-violation in B-meson decays is the principal goal of experiments at the current 
electron-positron B-factories, KEK-B and SLAC-B, and at the hadron colliders, Tevatron and LHC. In 
the standard model (SM), the source of CP viola tion is the Kobayashi-Maskawa phase^which resides in 
the Cabibbo-Kobayashi-Maskawa (CKM) matrixQEI. In the Wolfenstein parameterization^ of the CKM 
matrix, characterized by the parameters A, A, p and 77, CP violation is related to a non-zero value of 
the parameter 77. Of particular importance in the analysis of CP violation in the B-meson sector is the 
following unitarity relation: 

vm + VcdV*b + v td v t t = 0, (i) 

which is a triangle relation in the complex p — 77 plane, depicted in Fig. ^ Here, p = (1 — A 2 /2) p and 
fj = (1 — A 2 /2)?7 are the perturbatively improved Wolfenstein parameters^. The sides of this triangle, 
called Rb and R t , are defined as 



R b ee Rt = V(l-P) 2 +V 2 , (2) 

and its three inner angles have their usual definitions: 

QSMg HfeJ' 0s ^[-v^)- ' lsa ' g {-v^)- (3) 

The BELLE collaboration uses a different notation for these angles: cf>i = (3, cf>2 — a, and cf>3 — 7. We 
recall that among the CKM matrix elements above V u b and Vtd have sizable imaginary parts, and hence 
all three angles a, (3 and 7 are sizable. Of these, the phase (3 has already been well measured using the 
time-dependent CP asymmetries in the B — ► J/ipK~s and related decays, yielding El 

sin(2/3) = sin(20i) = 0.736 ±0.049, /3 = fa = 23.8° ± 2.0°. (4) 

The current thrust 3 of the two _B-factory experiments - BABAR and BELLE - is now on the measure- 
ments of the other two angles a (or §2) and 7 (or ^3). Of these, the weak phase a will be measured 
through the CP violation in the B — ► ttt:, B — > pn and B — > pp decays. To eliminate the hadronic uncer- 
tainties in the determination of a, an isospin analysis of these final states (as well as an angular analysis 
in the pp case) will be necessary^. To carry out the isospin analysis in B — > hit decays, one needs to know 
the three amplitudes A H , A 00 and A +0 , corresponding to the B® — ► ir + ir~ , B® — > 7r°7r° and B + — > 7r + 7r° 
decays, respectively, and their charge conjugates A v . At present, the only missing pieces in the current 
data are A 00 and A 00 - the amplitudes of the B® — > 7r°7r° and B^ — > 7r°7r° decays, respectively - though 
the measured charge conjugate averaged branching ratio ^ B(B°/B° — > 7r°7r°) provides an information 
on the sum |j4 00 | 2 + |j4 00 | 2 . Hence, a model-independent isospin analysis of the B — ► hit decays can not 
be carried out at present from the branching ratios alone. 

In addition to the measurements of the branching ratios in the B — > hit decays, the time-dependent 
CP asymmetries in the B®(t) / B®(t) — > ir + ir~ and B®(t)/B®(t) — > ir°ir° decays will greatly help in pinning 
down the weak phase a. We shall concentrate here on the CP-asymmetry in the decay B® — » ir + ir~, 
which is defined as follows: 

f ,s _ r[B° d (t)^n+n-]-r[B°(t)^« + «-] ^ 



(I 



+- 



r[BS(t)-f7r+7r-]+r[BS(t)-f7r+7r-] 

= S+- sm(AM B t) - C+- cos(AAf s t), 

where AMg is the mass difference in the B^ — B® system which is already well measured El and C£~ 
and S^T are the direct and mixing-induced CP asymmetry parameters, respectively. In the notation 
used by the BELLE collaboration WW C+~ is replaced by where A+~ = —C%~. 

The BABAREEl and BELLE WU measurements were summarized last summer at the Lepton-Photon 
2003 conference, yielding the world averages^ S+~ = —0.58 ±0.20 and C+~ = —0.38 ±0.16. However, 
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Figure 1 . The unitarity triangle with the unit base in the p — r\ plane. The two sides Rf, and Rt and the angles a (</>2 ) , /3 
(<j>l) and 7 ($3) are defined in Eqs. and respectively. 



a significant disagreement between the two measurements existed and the confidence level that the two 
are compatible with each other, in particular in the measurement of , was rather low (4.7% C.L.). 
Recently, the BELLE collaboration have updated their results for and A+~ by including more data. 
The current BELLE measurements ^21 (based on 140 fb _1 data) together with the updated BABAR 
results ® (based on 113 fb _1 data) of these quantities are as follows: 

„+__ J -0.40 ±0.22 ±0.03 (BABAR) +_ ( -0.19 ± 0.19 ± 0.05 (BABAR) , 
^ ~ \ -1.00 ±0.21 ±0.07 (BELLE) ' ™ ~ \ -0.58 ± 0.15 ± 0.07 (BELLE) ' ^ ' 

They have been averaged by the Heavy Flavor Averaging Group [HFAG] to yield ^1 

S+- = -0.74 ± 0.16, C+- = -0.46 ± 0.13, (7) 

and correspond to 4.6 and 3.5 standard deviation measurements from null results, respectively. It is also 
reassuring to note that the BELLE and BABAR measurements are now closer to each other than was the 
case at the Lepton-Photon 2003 conference, having now a scale factor e of 1.7 in S£~ and 1.4 in . 
Significant updates of the BABAR and BELLE results in the B — > tttt decays are awaited later this year 
which will further firm up these measurements. 

As can be judged from the results in (JJJ, current measurements of and Cj~ have already reached 
a significant level and invite a theoretical analysis leading to a determination of the unitarity triangle 
angles a and hence also 7. The importance of these measurements for 1 lie ('KM p henomenology has 
been long anticipated and discussed at great length in the literature IMLLMislLiial < Ug l -iy l Our analysis 
taking into account the updated B — > tttt data has many features which it shares conceptually with the 
cited literature and we shall compare our results with the ones obtained in the more recent works EISI 
A prerequisite to carry out such an analysis is to get model-independent bounds on the non-perturbative 
dynamical quantities \P C /T C \ and S c = 8^ — 6^, involving the so-called QCD-penguin P c and color- 
allowed tree T c topologies. Here, the subscripts denote that we are using the c-convention of Gronau and 
RosnerEHl in choosing the independent CKM factors in the analysis of the B — » tttt d ecays. Discussions of 
the ambiguities in the penguin amplitudes have also been presented Our approach makes 

use of the isospin-based bounds on the ratio \P C /T C \ and 5 C in the analysis of the data in the B — > tttt 
sector and we show how to incorporate these bounds in the analysis of the unitarity triangle in the SM. 

There are essentially three parameters \P C /T C \, 5 C and a [the weak phase (3 is already well measured, 
see Eq. @], which can not be determined from the measurements of just two quantities and C+~ . 
However, correlations and bounds on these parameters can be obtained which have been presented by 



e We thank the Heavy Flavor Averaging Group and, in particular, Andreas Hocker, for providing us the updated averages 
and the scale factors. 
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the BELLE collaboration based on their data EE! In the first part of our paper we undertake a similar 
analysis of the combined BABAR and BELLE data and work out the best-fit values and bounds on 
the parameters S c and \P C /T C \. As our analysis is performed within the SM, we allow the phase j3 to 
vary in the experimental range and restrict the range of a from the indirect unitarity-triangle (UT) 
analysis, which we have taken from the CKM fitterl=^and another recent fit of the CKM parameters^! 
We first show in this paper that the current data on S%~ and restricts the two strong interaction 
parameters 6 C and \P C /T C \. This information is already helpful in providing some discrimination on various 
competing approaches incorporating QCD dynamics in these decays. Conversely, restricting the allowed 
range of \P C /T C \ from the current dynamical models, data on S£~ and allows to put constraints 

on a. This has been done by the BELLE collaboration El, yielding at 95.5% C.L. 90° < cj> 2 < 146° for 
0.15 < \P C /T C \ < 0.45 and sin^^!) = 0.746. However, due to the restrictions on \P C /T C \, this remains a 
model-dependent enterprise. 

Our analysis differs in this respect from the one carried out by the BELLE collaboration. Instead of 
restricting \P C /T C \ by a survey of models, we use the isospin symmetry to restrict the range of \P C /T C \ 
and S c . To do this, we harness all the current data available on the branching ratios for B® — > ir + ir~, 
B + — > 7r + 7r and — > 7r°7r° decays (and their charge conjugates), S£~ and C+~ , and study a number of 
correlations, in particular |P C /T C | vs. cos(2#), where 9 is a penguin-related angle which is connected with 
the relative phase between the amplitudes and (see Fig. EJ. It is well known that the iso spin 
symmetry can be used to put a lower bound on cos (29), as first point ed o ut by Grossman and Quinn^l 
Subsequently, the Grossman-Quinn bound was improved by Charles 123 who derived in addition a new 
bound involving the i?° — > 7r°7r° and — > tt + tt~ decay modes. Based on the observation that the 
B — > 7T7T amplitudes can be rep rese nted, using the isospin symmetry, as two (closed) triangles which have 
a common base, Gronau et. al EH derived an improved lower bound on cos(2#) - the Gronau-London- 
Sinha-Sinha (GLSS) bound. We illustrate this bound numerically using current data and the constraints 
that it implies in the \P C /T C \ — cos(26') plane for both the 9 > and 9 < cases, varying 7 in a large 
range 25° < 7 < 75°, which adequately covers the present range of this angle allowed by the UT fits at 
95% C.L. Lest it be misunderstood, we emphasize that our final results for the CKM parameters and 
the dynamical quantities make no restrictions on the range of 7, whose value will be returned together 
with those of the other quantities by our CKM unitarity fits. The isospin-based lower bound on cos(2#), 
and hence an upper bound on \9\, is a model-independent constraint on the penguin contribution in the 
analysis of the data involving the measurements of S£~ and C+,T . 

There are yet other bounds based on the isospin symmetry in B — > 7T7t decays which lead to restrictions 
on 7. In parti cular, the Buchalla-Safir bound "™ on 7 (and its various reincarnations discussed recently in 
the literature EHBj ) result from the correlations involving sin(2/3), 7, S^T, and C+,7. We have analyzed 
these bounds, but we find that they are not very useful at present as the current central values of sin(2/3) 
and — S+~ almost coincide. For these bounds to be useful phenomenologically, the value of —S+~ has to 
come down substantially. 

In the last part of our analysis, we study the impact of the and C+^T measurements on the profile 
of the unitarity triangle in a model-independent way. We first show that the quality of the UT fits is not 
modified by the inclusion of the data on and C£~ , as the two additional parameters \P C /T C \ and 5 C , 
when varied in large regions, can always reproduce the central values of the and averages. 

We then implement the lower bound on cos(2#) in performing the fits of the unitarity triangle in the 
p — fj plane. The present bound cos(2#) > 0.27 removes a small part of the otherwise allowed region 
of the unitarity triangle, but this constraint will become more significant in future as the errors on the 
B — > 7T7T branching ratios, S+~ and are reduced. The effects of the bound on cos(20) are also shown 
on the correlations a — 7 and cos(2a) — cos(2/3). Working out the ^-distributions in the quantities a, 
\P C /T C \ and S c , we find that the current data prefers rather large values for the latter two quantities, 
with the minimum of the ^-distributions being at \P C /T C \ = 0.77 and S c = —43°. The corresponding 
best-fit values of a and 7 are a = 92° and 7 = 64°. At 68% C.L., the ranges are: 81° < a < 103°, 
21.9° < (3 < 25.5°, 54° < 7 < 75°, 0.43 < \P C /T C \ < 1.35, and -64° < S c < -29°. The bound on cos(26>) 
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is very efficient in the exclusion of the large values of \P C /T C \ and —6 C . Their best-fit values are quite 
a bit larger than anticipated in most dynamical approaches. This feature has also been noted in earlier 
studies on the B — > -kit decays E3I±E§I 

This paper is organized as follows: In section 2, we give the relations among the observables in 
the B® - * 7r + 7r~ decay and its charge conjugate, the CKM parameters and various dynamical quantities. 
Section 3 contains a review of several isospin-based bounds in the B — > 7T7t decays. In section 4, we report 
on the results of our numerical analysis of the time-dependent CP asymmetry in the B^/B^ — ► 7r + 7r~ 
decays, and in section 5 we show the results of the unitarity triangle fits, correlations involving the 
angles a, f3 and 7, and the dynamical quantities \P C /T C \ and 5 C , carried out in the context of the SM. We 
conclude with a summary and some remarks in section 6. 



2 Relations among the Observables in the B — ► tttt Decays, CKM Parameters and 
Dynamical Quantities 

In this section we present the analytic formulae that we need to discuss the time-dependent CP asymmetry 
and the branching ratio in the — > ir + ir~ decay, and their relations with the CKM parameters and 
various dynamical quantities. 

The amplitudes = A[B^ — > tt + tt~] and its charge-conjugate A^ = A[B^ — > ir + ir~] can be 

written by using the Gronau-Rosncr c-conventionl^ as follows: 

a + - = v; b v ud a+- + v; b v cd a+- + v; b v td a+- (8) 

= V* b V ud (A+- - A+-) + V* b V cd (Af- - A+-) 
= -(|T c |e< e +<7 + |p c | e *f ) , 
A+- = V ub v: d {A+- A+-) + V cb V c * d (A^ - Af~) (9) 
= - (|T C | e i5 ° e-^ + \P C \ e<) . 

In getting the last expressions for the amplitudes the unitarity relation has been used together with 
the phase convention V ub = \V u b\ e~' n for this CKM matrix element. 

The phenomenon of the B® — B® mixing modulates the time dependence of the decay amplitudes for 
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B d(t) 7r+7r_ and 7r+7r_ 



= e -iM B t e -rt/2 



AM B t 



cosh — ■ Ai_ sinh — ■ — 



I sin ■ 



AM B t 



AlL cosh — — - — sinh — — 



- e -iM B t e -Tt/2 



COS ■ 



4 

AM B t 



.4 4 



A+ cosh ^— ^ — sinh ^— ^ 



(10) 



(11) 



z sin ■ 



AM B t 



cosh ^— ^ — Air sinh 



P -A+- 



Here, M B and T are the average mass and decay width of the — B^ system, and AM B and Ar are 
the the mass- and width- difference in the two mass eigenstates, respectively, p/q ~ V t * d /Vtd = e 2 ^ is the 
mixing parameter, and the quantity: 



q A+- 



-Ti 



-Tt 



\+- - 1 - P 2i« z 1 c/ c| = |a+-| e 2iQrff (12) 

P A+- 1 + \P C /T C \ e iS ° e-*t ~ 1 ™ 1 ' { ' 

is introduced which encodes all the information about the CP asymmetry in this decay. Here, a e s = a + 9, 
and 9 is the penguin-pollution parameter shown in Fig. which is connected with the relative phase 
between the amplitudes and A A , A</> H = 2(7 + 6*), and the relation a + /3 + j = tt has been used. 
Note that in the limit P c /T c — > 0, 9 — > and a c g — ► a. 

The partial decay widths of the time-dependent B%(t) — > 7r + 7r~ and — * 7r + 7r~ decays are 
proportional, respectively, to^^ 

(13) 
(14) 

(15) 
(16) 

(17) 

(18) 



\A+-(t)f 
\A+-(t)f 



B+~ [1 + C+~ cos(AAf B i) - S+- sin(AAf B <)] 
B+- [1 - G+T cos(AAf B i) + 5+T sin(AAf B <)] 



where |p/?| = 1 is used and the following quantities are introduced: 



b; 



\A 



+ -|2 



\A+- 



= |A+-| 2 -|A+-| 2 = 1-|A+-| 2 
™ |A+-| 2 + |I+-p 1 + |A+7| 2 ' 



^tir — 



Vtir 



2lm[(q/p)A+-(A+-y] 2ImA+- + 

\a+^TWT 2 = TTp4fF = y ™ ( eff) ' 

2|A+-||A+-| 2|A+"| 



|4+-| 2 + |A+-| 2 1 + IA+7I 2 ' 

Using the expression for given in Eq. I|12|) . the above quantities can be rewritten in the following 

form: 



|T c | 2 i?+- - |T c | 2 +2|P c ||r c |cos5 c cos 7 +|P c | 2 , 



Si 



+- 

7T7T 



Rtttt 

1 



sin 5 C sin 7, 



sin(2a) - 2 



cos 5 C sin(a — (3) 



sin(2/3) 



Making the back transformation, \T C \, \P C \ and S c can be expressed as follows: 



Bt£ 
1 - cos(27) 



[l-y+-cos(20-2 7 )] 



(19) 
(20) 

(21) 

(22) 

(23) 



l^| 2 = IZ^)[l-^-cos(2.)], (24) 

tan^ c = + ^ Si 7 • (25) 

yj n cos(26> - 7) - cos 7 

In the limit of neglecting the penguin contribution (i.e., \P C \ — * 0), the CP-asymmetry coemcient C+~ 
goes to zero (and — > 1) as well as 9 — > 0, in agreement with Eq. I|24|) . Also, in this limit, S+^r = |T C | 2 , 
as in this case the branching ratio is completely defined by the tree contribution. 

In terms of 7 and 9, the penguin-to-tree ratio squared has the following expression: 

2 

„2 



1 ~ vtn cos(26>) = 1 - yj v cos(2o; cff - 2a) ^ 



1 - yt-x cos(26» - 27) 1 - y£„ cos(2a e ff + 2/3) 

This relation constrains r c in terms of cos(2#), given 7 and y^~ ' . It should be noted that for fixed values 
of y+~ and r c , cos(27) varies in the range: 

1 - ( C+-) 2 (l + ^ )/(2rg) 

1 (CtJtt )^ 

It is easy to see that the upper limit of cos(27) is equal to 1 when r c — 1, independent of . For 
r c 7^ 1, the allowed range of cos(27) puts a constraint in the r c — plane. Thus, the allowed domain 
of r c is completely defined by the direct CP-asymmetry coefficient and, for negative (in accordance 
with the experimental data), it is given by 



1 < cos(2 7 ) < K ™ ' cJ . (27) 



i±^i-(c+r) 2 

(^c)min.max — ■ (28) 

In particular, for the central experimental value = —0.46, the allowed range of r c is as follows: 

0.244 < r c < 4.104. (29) 

For the current experimental central value of C£~ , and its ±1<7 limits, the dependence of the upper limit 
cos(27)|ul on the magnitude of the penguin-to-tree ratio is presented in Fig. El Decreasing the magnitude 
of , the allowed region of r c becomes wider. The expression for cos(26*) as a function of y^T, r c and 
cos(27) is as follows: 



(1 - r c 2 )[l - r c 2 cos(2 7 )] ± r c 2 J [I - cos 2 (2 7 )]{2(y+ 7r -) 2 r 2 [l - cos(2 7 )] - (1 - r c 2 ) 2 [l - (y+T) 2 ]} 

cos 26» = ¥— r . 

1 ' y+r{(l-r 2 ) 2 + 2r 2 [l-cos(2 7 )]} 

(30) 

3 Isospin-Based Bounds in B — > tttt Decays 

With partial experimental information on B — > tttt decays available at present, it is of practical importance 
to get useful restrictions on the dynamical parameters at hand, r c and 5 C . It is obvious from H25|) and l|26|). 
that apart from Cj~ , which is measured from the time-dependent CP asymmetry, the angle 29 plays a 
central role in constraining the dynamical parameters of interest. While 29 will be determined eventually 
from the measurement of A 00 and A 00 (see Fig.|2J), this is not the case now. Instead, several bounds have 
been derived on cos(2#) using the isospin sy mme try, which we review here. The first of these which we 
will work out numerically is the GLSS bound ^21 

( b+- + 2B +0 - 2B Q0 ) 2 - 4B+-B+ 

COs(2fl) > t z Tn ) (oil 

^yn-n j-'tttt ^ittt 

where B+® and B®^ are the quantities constructed from the B + — ► tt + tt° and B® — ► 7r°7r° decay ampli- 
tudes in a similar way as in Eq. I|15[) . It was also demonstrated by Gronau et. al ESI that this bound is 



7 




Figure 3. The dependence of the upper limit cos(27)tjl on the ratio \P C /T C \ for Cjt n = —0.46 (the current central value), 
= —0.59 and C^tT = —0.33, which demarcate the ±ler experimental measurements. 



stronger than both the Grossman-Quinn and Charles bounds. As a byproduct, a bound on the direct CP 
asymmetry in the — > 7r°7r° decay was also obtained by these authors ESI 



B + ~ ( B+- - 2 R+° - 2 R 00 "1 



<^fUU \ /-^fH 7T7T \ 7T7r 7T7T TTTT J ('^9) 

" ^ 2i?™ (S+T + 2B+° - 2B°° ) • 1 j 

There have also been attempts to derive isospin bounds on cos(2(9) in the £>° — > 7r + 7r~ decay 
which are based on the knowledge of the direct CP asymmetry alone. As for the GLSS bound ^"1 
the starting point is the isospin-based triangular relation between the , A 00 and A +0 amplitudes: 

A+- 



V2 



A 00 = A +0 , (33) 



corresponding to the — > 7r + 7r~, B^ — > 7r°7r° and B + — ► 7r + 7r° decays, respectively, and a similar 
one for the phase-shifted charged-conjugate amplitudes A 1 ^ = e 2 * 7 ^4 l - J . The graphical representation 
of both triangles is shown in Fig. [21 It should be noted that the magnitude of the difference between 
the A^ and A^ amplitudes is -s/2|i-c| sin 7 and not \/2| -Ft | sin a as the c-conventionGOl is employed for 
the amplitudes throughout this paper. With the help of the sine theorem, sin \29\ can be written as: 

sin|20| = ^psin^ = «^ S in^. (34) 



Squaring all the terms, the above relation can be rewritten as two inequalities: 

S in 2 l2g]<2 1 -^ C ° S _ (2 ^ ) s^\2e\<2 l - y ™ CO S 6 \ (35) 



following from the conditions sin 2 9 a < 1 and sin 2 9a < 1, respectively. Here, Eqs. (|15|l . Ijltil) and (|24|l 
were used to eliminate \A^ | 2 , \A^ | 2 and |P C | 2 . While sin 2 \29\ < 1, this does not imply that the 
expressions on the r.h.s. of l|35|l also satisfy this upper bound. Hence, no bound on \29\ follows from 
Eq. Qgftf and the GLSS bounds are indeed the strongest isospin-based bounds in the B — > tttt sector. 



We are grateful to David London and Nita and Rahul Sinha for pointing this out to us. 
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In addition to the above bounds on cos(2#), bounds on the CKM angle 7 have also been derived in 
the literature recently which are based on the study of the correlation 7 — , given sin(2/3) | 3U | 31 | 19 | 
We reproduce these bounds below and discuss their impact in the next section. Relating the unitarity 
triangle angles (3 and 7 with the Wolfenstein parameters p and fj: 

1- p±ifj = R t e ±l(3 , p±ifj = R b e ±l1 , (36) 

where R t and Rb are defined in Eq. @ , the quantities , C£~ and 5 + „r can be expressed in the form: 



'■7T7T 



' TV 



Rb 




T c 


2fj 




RbRtn 




-2fj 




p2 d2 pH — 



cos S c + 
sin <L , 



p-Rl + {l-Rl)R b 



COS<5 c + (1 - ,o)i? 2 





2" 


T c 





(37) 
(38) 

(39) 



The relation for given above agrees with Eq. (5) of the paper by Buchalla and Safir Is"*! if 

introduces the pure strong-interaction quantity 



r = R b 



Pc 

T, 



(40) 



used by these authors. Note also that the equation for can be rewritten in terms of this quantity r 
in the following form: 

2 / ^_ x 2 



(p + r cos<5 c ) + rj 



r sin 5 C 

r+- 



r sin S, 



(41) 



For the phenomenological analysis, it is more convenient to eliminate p from Eqs. (|38(l and Q39JI with the 
help of the relational 



1 — p — r\ cot /3 = r)T. 

With this, the Wolfenstein parameter fj can be related to either S£~ or C£~ as follows: 

1 



j) 



(S) 



(l + r 2 )5+.- 



(l + TS+-)(l + rcos«5 c ) 



± yj(l - 5+T) 2 (l + r cos<5 c ) 2 - (1 + t 2 )5+7 [S+T + sin(2/3)]r 2 sin 2 <J C 

1 



n {C) = J(l + rcos5 c )sin(2/3) + - , 

i (l + T 4 jG,r, 



• sin <5 r 



(42) 



(43) 



(44) 



± y (1 + T 2 )(yt^ ) 2 r 2 sin 2 5 C - [C^ (1 + r cosd c ) - Trsm6 c ]' 

The first of these relations has been obtained by Buchalla and Safir (BS) 2D| and has been used to derive 
an upper bound on fj and, hence, a lower bound on 7: 



7 > arctan ■ 

' ~ 2 



1- \/i-(s+r) 2 



1 + tS+T-J1-(S. 



(45) 



which holds in the range — sin(2/3) < Sj w < 1. However, the current central experimental value J7J 
of practically coincides with the central value Q of — sin(2/3), and hence no useful bound on the 



CKM angle 7 follows from the BS bound at present. We shall show this bound as a function of S+ v , as 
well as its extension for the case of 7^ 0: 



1 + S+- sin(2/3) + Jl - (C+T) 2 - {S+«Y cos(2/3) 

tan 7 >L_ = v , (46) 

yj 1 - (C^-) 2 - (S+„) 2 sin(2/3) - 5+T cos(2/J) 

obtained by Botella and Silva l&M The next step in the generalization of the BS bound was recently 
undertaken by Lavoura ^ who considered the modification of this bound by putting restrictions on the 
strong phase 5 C . With the current experimental values of and C+~ , also the Botella-Silva and 
Lavoura versions of the BS bound are currently not useful in constraining 7. With precise measurements 
of and in future, these bounds may, in any case, provide useful consistency checks for the 

dynamical models used in the estimates of r and S c . 



4 Numerical Analysis of the CP Asymmetry in — * tt + tt Decay 

Within the SM, the targets for the experiments measuring the angles a and 7 are fairly well defined, 
as the fits of the unitarity triangle through the measurements of the CKM matrix elements yield the 
following ranges for these angles at 95% C.L.E3 

70° < a < 115° , 43° < 7 < 86° . (47) 

The corresponding 95% C.L. ranges obtained using the default values of the input parameters by the 
CKM fitter group are very similar: 

77° < a < 122° , 37° < 7 < 80° . (48) 

So, if the SM is correct, and currently there is no experimental reason to believe otherwise, then from the 
B — > 7T7T analysis, values of a and 7 should emerge which are compatible with their anticipated ranges 
listed above. Of course, the hope is that direct measurements of these angles will greatly reduce the cur- 
rently allowed ranges. However, for this to happen, one has to determine the dynamical quantities \P C /T C \ 
and 5 C . 

Surveying the recent literature on the estimates of \P C /T C \ and 5 C in B — > 7T7t decays, we remark 
that they are either based on specific schemes based on facto rizatio n in which non-factorizing effects are 
implemented using perturbative QCD in the large-m?, limit l ^ l ^ l or on phenomenological approaches 
based on some input from other data and factorization. A typical study in the latter case makes use 
of the data on the B — > n£vg and B — > Kit decays, which a re u sed in conjunction with the assumption 
of factorization and estimates of the 5J7(3)-breaking effects Some representative estimates in these 
approaches are as follows: \P C /T C \ = 0.285 ± 0.076 [Beneke, Buchalla, Neubert, Sachrajda] ESI \P C /T C \ = 
0.32tg;i£ [Beneke, Neubert] El, \P C /T C \ = 0.29 ± 0.09 [Buchalla, Safir] 30 , \P C /T C \ = 0.23±S;g£ [Keum, 
Sanda]EBl \P,/T C \ = 0.276 ±0.064 [Gronau, Rosner] 7 , \P C /T C \ = 0.26 ±0.08 [Luo, Rosner]E3. (See, also 
Xiao et al. I***). Thus, \P C /T C \ = 0.30 is a typical value from these estimates. 

What concerns the strong phase difference S c , the two dynamical approaches developed in detail 
(QCD-Factorization and pQCD El) differ considerably from each other due to a different power 
counting and the treatment of the annihilation contributions in the decay amplitudes. When comparing 
the current data with these specific approaches, we shall take for the sake of definitencss the estimates 
by Buchalla and Safir ESI to represent the QCD-factorization approach, \P C /T C \ = 0.29 ± 0.09 and 6 C = 
0.15 ± 0.25 radians (S c = 9° ± 15°), and the estimates by Keum and SandaESI \P C /T C \ = 0.23^5 and 
—41° < S c < —32°, for the pQCD approach. Within the SM, the consistency test of these approaches 
lies in an adequate description of the data on and C+^T, with the parameters a, 7, \P C /T C \ and S c 
all lying in their specified ranges. However, as |P C /T C | and 5 C are not known directly from data or a first 
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principle calculation, we can leave them as free parameters and determine them from the overall fits. We 
shall pursue both approaches in this section. 

We now present our numerical analysis of the current averages of and given in (0. For 

the construction of the C.L. contours, the following x 2 -function is used: 



x 2 = 



S n7r (^S 7rn ) 



cxp 



(49) 



which is equated to 2.30, 6.18, and 11.83, corresponding to 68.3%, 95.5%, and 99.7% C.L., respectively, 
for two degrees of freedom. 

We start by showing that the current data on and in the B® — > tt + tt~ decays provides a 
discrimination among various dynamical approaches, for which the QCD factorization ^ and perturbative 
QCD ° approaches will be taken as the two leading contenders. The results of this analysis are presented 
in Fig. 0] for six values of a in the range 80° < a < 130° in intervals of 10°. To take into account the 
dispersion in the values of \P C /T C \, we take three values of this ratio, namely 0.30, 0.55, and 0.80. The 
first of these values represents the current expectations of this quantity, whereas the last is taken with 
the hindsight of the best fit of the data that we have performed in a model-independent way, as described 
later. The points indicated on these contours represent the values of the strong phase difference S c which 
is varied in the interval — tt < S c < tt. We do not show the plot for a = 70°, which is the 95% C.L. lower 
value of a from the unitarity fits, as already the case a = 80° requires rather large value of \P C /T C \. In 
each figure, the outer circle corresponds to the constraint (S+~) 2 + (C^r) 2 = 1. The current average JZj) 
of the BABAR and BELLE data satisfies this constraint as shown by the data point with (unsealed) 
errors. The two ellipses surrounding the experimental measurement represent the 68.3% and 95.5% C.L. 
contours. This figure demonstrates that, as C+~ is negative and large, current data favors a rather large 
strong phase, typically —60° < S c < —30°. The two shaded regions shown in this figure correspond to 
the predictions of the QCD-factorization approach (the upper shaded area) and the perturbative QCD 
framework (the lower shaded area). As can be seen, the predictions of the QCD-factorization approach lie 
outside of the 3a exp erimental measurements for all values of a shown in this figure. For the perturbative 
QCD framework^!, one finds agreement with the measurements, but only at about 2a level. Restricting 
a in the region 90° < a < 110°, good fits of the data are obtained for typically \P C /T C \ > 0.5 and 
5 C < —30°. We shall quantify the fits more precisely later. 

We now turn to a model-independent analysis of the and data. As the current world aver- 
ages of these quantities are negative and rather large (7J, positive values of the strong phase difference S c 
are excluded at a high confidence level (> 99.7% C.L.) which is demonstrated in Figs. 03 Eland[7| This 
observation is the reason why the full range — tt < 5 C < +tt has been restricted to the negative values of S c , 
— 7r < 5 C < 0, in these figures. In Fig. [SJ three representative values (90°, 105° and 120°) of the angle a 
from the UT-favored interval l|48|l and three values (25.9°, 23.8° and 21.9°) of the angle f3, which cover 
the present measurement of this quantity within ±ler range are shown and the resulting x 2 -contours 
in the variables \P C /T C \ and S c are plotted. Note that the dependence on the precise value of (3 in the cur- 
rent experimental range of this angle is rather weak. Hence, we show the /3-dependence of the correlation 
for only one value of a, namely a — 105°. The most important message from this analysis is that the 
current data favours negative and rather large values of the strong phase 5 C , which are correlated with 
the values of a. Restricting to the 68.3% C.L. contours for the sake of definiteness, the minimum allowed 
values of — S c are: 30°, 45°, and 70° for a = 90°, 105°, and 120°, respectively. What concerns the allowed 
values of \P C /T C \, we note that except for a relatively small allowed region near 80° < a < 90°, they 
overlap with the theoretical estimates of the same specified above at 95.5% C.L. However, the best-fit 
values of \P C /T C \ are on the higher side as shown by the dots in these figures. It should be noted that 
the current data results in the lower bound on the penguin-to-tree ratio \P C /T C \ > 0.18 at 95.5% C.L but 
extends to much larger values of \P C /T C \, which we have suppressed in these figures for the sake of clarity 
but will show in the next section where we discuss the fits of the unitarity triangle. 

The correlations between a and S c for three fixed values \P C /T C \ = 0.25, 0.35, and 0.45 in the 
theoretically motivated interval are shown in Fig. This figure updates the results by the BELLE 
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Figure 4. Implications of the time-dependent CP asymmetry parameters C^tT = —0.46 ± 0.13 and = —0.74 ± 0.16 

from the BELLE and BABAR, measurements for the CP violating phase a (or <f>2). In this analysis, the strong phase 8 C is 
varied over the full range —it < 8 c < tt and curves are drawn for three values \P C /T C \ = 0.30, 0.55, and 0.80. The predictions 
of the QCD factorization (upper box) and pQCD (lower box) approaches are also shown for fixed values of a noted on the 
six frames. The curves around the data point represent the 68.3% and 95.5% C.L. contours. 



12 



collaboration and shows that a satisfactory description of the current data for these values of \P C /T C \ 
and with a lying within the indirect UT-based range is possible only with large values of the strong 
phase —5 C . Again using the 68.3% C.L. contours, it is seen that the minimum allowed value for |P C /T C | = 
0.35 is -8 C ~ 55°, and it decreases to 45° for \P C /T C \ = 0.45. With 0.20 < \P C /T C \ < 0.45, the angles a 
and 8 C lie in the intervals: 90° < a < 130° and -160° < 8 C < -30°, at the 95.5% C.L. As higher values 
of \P C /T C \ are experimentally allowed, the correlations between a and 8 C for larger values of \P C /T C \ = 0.55, 
0.65, and 0.75 are shown in Fig.0 We note that with these values, the allowed ranges for the angles a 
and 5 C become wider with increasing \P C /T C \. 

The correlations between the angle a and \P C /T C \, for three representative values of the strong-phase 
difference 6 C = —40°, —80°, and —120° and the angle (3 within its experimental range, are presented in 
Fig- ID This figure demonstrates again that smaller values of |5 C | require larger values of \P C /T C \. The 
restrictions on \P C /T C \ and a discussed above are also seen in this figure. 

In summary, we see that current data allows a wide range of the quantities \P C /T C \ and S c , and 
without restricting them the impact of the and measurements on the CKM parameters, in 

particular the angles a or 7, is rather small. The dynamical approaches discussed above are not a great 
help as they are not good fits of the data within the SM. 



4--1 Constraints on cos(2#) from bounds based on the isospin symmetry 

As discussed in the previous section, the penguin contribution in the B® — ► ir + ir~ decay can be param- 
eterized by the angle 29. Having at hand the experimental range of , it is of interest to work out 
numerically the dependence between the ratio \P C /T C \ and cos(2#), given by Eq. Ij26(l . The results of this 
analysis are presented in Fig. [5] for six values of 7 in the range 25° < a < 75° in intervals of 10°. The 
solid lines in all the frames correspond to the central experimental value of while the dashed 

lines correspond to the ±1<7 values of this quantity. Due to the functional dependence (|26|l of r\ on 
cos(2# — 27), there exists a sign ambiguity and there are two solutions depending on 9 > and 9 < 0. We 
show both of these solutions and each frame in this figure contains two sets of curves where the upper 
and the lower ones correspond to 9 > and 9 < 0, respectively. The isospin symmetry and the existing 
data on the B — > irir decays allow to put restrictions on cos(2(9). The GLSS lower bound H31|l on cos(26>) 
is based on the B — > ixix branching ratios and y~^~ . The recent experimental data on the branching ratios 
and the B + - and £>°-meson lifetime ratio ^3 

B{B° d -> tt+tT) = (4.55 ± 0.44) x 10~ 6 , 
B(B + -» tt+tt ) = (5.27 ± 0.79) x 10~ 6 , 
B{B° -> 7r°7r°) = (1.90 ± 0.47) x 10~ 6 , 
r B +/r B « = 1.086 ±0.017, 

have been used in getting the conservative numerical bound: cos(2#) > —0.03. This bound is shown as 
vertical dashed lines in all the frames in Fig. [5J It should be noted that if the central values of the data 
are used instead, the resulting GLSS bound is: 

cos(20) > 0.27 , (50) 

GLSS 

which is shown as the solid vertical lines in Fig. |5| The shift is mainly due to the current uncertainties 
in the branching ratios for B® — > 7r 7r° and B + — > 7r + 7r°. Our analysis shows that putting a lower bound 
on cos(2#), \P C /T C \ gets significantly constrained. It is seen that the branch with 9 < results in smaller 
values for \P C /T C \, which are concentrated in a relatively narrow interval. However, as 7 increases, this 
interval becomes wider. A priori, it is difficult to argue which of the two solutions 9 > and 9 < should 
be entertained. Hence, in the implementation of the isospin-based bound on cos(2#) in the unitarity fits, 
we shall allow the sign of 9 to take either value. 
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Figure 5. The correlation \P C /T C \ — 8 C corresponding to the 68.3%, 95.5%, and 99.7% C.L. ranges of C^n and S^n for 
three values of a and /3. Note that the dependence on /3 is rather weak and hence not shown for the other two values of a. 



Based on the central values of the experimental data specified above and C+ v Q , the minimal value 
of the direct CP asymmetry in the — > ir°ir° decay (13211 can be estimated as: 



C°° > 0.47. 



(51) 



It should be noted that differs in sign from . (See, also the recent analysis by Buras et al. ^1) 
The SM-based bounds on the angle 7 as a function of the CP asymmetry with = (the 
Buchalla-Safir bound EH) and with C+~ = -0.46 ± 0.13 (the Botella-Silva bound E2) are shown in 
Fig. EI] as the solid line and the shaded area, respectively. The vertical band corresponds to the current 
experimentally measured value, and the central value practically coincides with S+~ — ~ sin(2/3). Thus, 
these limits do not provide any restrictions on 7 at present, but if with improved data a sizable shift of 
the central value from its current value takes place, then these bounds may lead to useful constraints. 
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5 Analysis of the CKM Unitarity Triangle Including and Measurements 



In this section we investigate the impact of the and measurements on the unitarity triangle 

fits. We adopt a bayesian analysis method to fit the data. Systematic and statistical errors are combined 
in quadrature. We add a contribution to the chi-square for each of the inputs pres ented in Table ^ 
Other input quantities are taken from their central values given in the PDG review^! The lower bound 
on AMb b is implemented using the modified-x 2 method (as described in the CERN CKM Workshop 
proceedings which makes use of the amplitude technique E3I The B s B s oscillation probabilities 
are modified to have the dependence P(B S — > B s ) oc [1 + _4cos(AAfs s t)} and P(B S — > B s ) oc [1 — 
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Acos(AMb b t)]. The contribution to the x 2 -function is then 

X 2 (AM B J = 2 



Erfc" 



- Erfc — = 

2 



-1 2 



(52) 



where A and a a are the world average amplitude and error, respectively. The measurements of 
and S^JT contribute to the x 2_ function according to Eq. (|49|l . The resulting x 2 -function is then minimized 
over the following parameters: p, fj, A, B K , rji, rj 2l m c (m c ), m t (m t ), r] B , /s d y/ B Bd , f, \P C /T C \, and 5 C . 
Further details can be found in Ref. ESI 

We present the output of the fits in Table where we show the 68% C.L. ranges for the CKM 
parameters, the angles of the unitarity triangle, AMb„, \P c /T c \ and S c . Note the enormous ranges 
for \P C /T C \ and S c allowed by the UT fits. The 95% C.L. constraints from the five individual quantities 
(Rb, £k, AMs d , AMb s , and a^K s ) an d the resulting fit region (the shaded area) are shown in Fig. 1111 
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Further details of this analysis and the discussion of the input parameters can be seen elsewhere^! The 
shaded areas in Fig.^]are the 95% C.L. correlations between a — 7 (the left frame) and sin(2/3) — sin(2a) 
(the right frame). In Fig. ^5] we show the behaviour of Xmin as a function of the angle a (the dashed 
curve). The solid curves in all these figures will be explained below. 

Note that the inclusion of the C~%~ and S£~ measurements does not induce any additional constraint 
on the fits. This is because we added two additional terms to the x 2 -fu nc tion with the dependence on 
two more variables \P C /T C \ and 5 C . Indeed, for any value of p and f), it is always possible to choose \P C /T C \ 
and S c so as to exactly reproduce the C+~ and S+~ experimental central values. Thus, the total % 2 of 
the unitarity triangle fit remains unchanged as the new measurements do not contribute to the total x 2 - 
In the two plots presented in Fig.^Jwe fix \P C /T C \ (the left frame) and 5 C (the right frame) and minimize 
the % 2 -function with respect to all the other variables. In both cases, the absolute minimum of the curve 
coincides with the minimum \ 2 °f the overall fit (Xmin ~ 0.57). A peculiar feature is the presence of two 



17 





0.0 



0.0 



-1.0 



-0.6 



-0.2 0.2 

cos(20) 



0.6 



1.0 



-1.0 



-0.0 



-0.2 0.2 

cos(20) 



0.0 



1.0 




Figure 9. The correlation |P C /T C | — cos(20) based on the current measurements of C^tT for fixed values of the angle 7 
indicated on the frames. The red dashed and green solid vertical lines represent the Gronau-London-Sinha-Sinha (GLSS) 
lower bounds cos(2(9) > —0.03 and cos(26>) > 0.27, respectively, as discussed in the text. The upper and lower sets of curves 
correspond to the branches with 8 > and 8 < 0, respectively. 
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Figure 10. The SM limits on the angle 7 in dependence on SiiT at C^T = (the Buchalla-Safir limit) and at CiT = 
—0.46 ±0.13 (the Botella-Silva limit) shown as the solid line and the shaded area, respectively. Note that both are the lower 
limit for Si^~ > — sin(2/3) and the upper one for S^T < — sin(2/3). The vertical band corresponds to the experimentally 
measured value of S^k with the recent central value which practically coincides with S^T = — sin(2/3). 

distinct regions for which the overall \ 2 is very small (the dashed curves). The best-fit values for \P C /T C \ 
and 5 C are \P C /T C \ — 0.77 and 5 C = —43°, respectively. Requiring higher confidence levels for the fits we 
obtain: 



Up to now we have not considered the impact of the isospin-based bounds in the analysis of the B — > 
7T7T data on the fits of the unitarity triangle. In future, the implementation of these bounds could take 
the form of the GLSS bound shown in Fig. [5] and the BS bound shown in Fig.^3 However, due to the 
proximity of and — sin(2/3) with each other in the current data, we will concentrate on implementing 
the GLSS lower bound on cos(2#) resulting from the isospin-based analysis of the B — > tttt data presented 
in the previous section. A possible implementation of this bound could be undertaken by minimizing the 
X 2 -function and rejecting points for which cos(2f?) lies below the allowed range. Taking cos(20) > 0.27, 
this analysis results in the black contours in Figs. II II and 1121 Note that a part of the 95% C.L. region is 
now excluded. The impact of the cos(2f?) lower bound on the x 2_n t f° r the angle a is shown through the 
solid curve in Fig. ED This does not change the best-fit value of a but restricts the allowed range of a by 
a few degrees at 95% C.L. The impact of the cos(2#) lower bound on the fits of the unitarity triangle and 
the correlations a — 7 and sin(2a) — sin(2/3) is currently not great, but this will change with improved 
measurements leading to tighter constraints on cos(2#), and eventually its measurement. The effect of 
the lower bound on cos(2#) is, however, very significant for the allowed value of \P C /T C \ and 5 C . This is 
shown in Fig. 1141 through the solid curves for \P C /T C \ (left frame) and S c (right frame). While the best-fit 
values of these parameters have not changed, the allowed regions are now drastically reduced. Thus, at 
68% C.L., the allowed values are: 



assumption. Our results in (|55|l can be compared with the analysis by Buras et al.<^^, obtained in the SM 



\P C /T C \ > 0.23 & S c < -18° @ 95% C.L. , 
\P C /T C \ > 0.15 & 5 C < -13° @ 99% C.L. . 



(53) 
(54) 
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by restricting (3 and 7 in the ranges 2/3 = (47 ± 4)° and 7 = (65 ± 7)°, which yields \P C /T C \ = 0.49±g;|? 
and S c — (—43123)°- The two analyses are compatible with each other though they differ in the details, in 
how the exact isospin-relations were imposed in the analysis of the data (and also somewhat in the input 
data). However, we note that we have not restricted 7 to any range, as it is a fit parameter returned by 
th e un itarity fits, but our fit value 7 = (64 ± 10)° is compatible with the input value used by Buras et 
alM 

Finally, to show the impact of the cos(20) bound on the unitarity triangle in a Gedanken experiment 
where the quantities S+~ and are assumed to be very precisely measured, we fix S£~ and C+~ to 
their current experimental central values JJJ and show the allowed region in Fig. 1151 resulting from the 
lower bound cos(26>) > 0.25 (the shaded region). Note that this results in a constraint in the p — fj plane 
which is very similar to what one gets using a range for a. Since the experiments do not measure a, but 
rather S£~ and Cj~ , this figure shows how the eventual S£~ and Cj^T measurements together with the 
lower bound on cos^S) gets translated. This represents the strongest possible constraint on the profile 
of the unitarity triangle from the and measurements that one can get in a model-independent 
way. Of course, this figure itself is only illustrative, as the actual constraints will depend on the values 
of S£~ and C£~ and the lower bound on cos(26>) that will be eventually measured. 



A 0.2224 ±0.002 (fixed) 

\V cb \ (41.2 ± 2.1) x 10- 3 

\V ub \ (3.90 ± 0.55) x 10" 3 

a^K s 0.736 ±0.049 

\e K \ (2.280 ± 0.13) x 10- 3 

AM Bd (0.503 ± 0.006) ps" 1 

7?i(m c (m c ) = 1.30 GeV) 1.32 ± 0.32 

t] 2 0.57 ±0.01 

rj 3 0.47 ±0.05 

m c (m c ) (1.25 ±0.10) GeV 

m t (m t ) (165 ±5) GeV 

B K 0.86 ±0.15 

f Bd ^B^ d (215 ±11 ±15123) MeV 

T] B 0.55 ±0.01 

£ 1.14 ± 0.03 ±0.02l°-J 3 l^ 3 

AM Bs > 14.4 ps- 1 at 95% C.L. 



Table 1. The input parameters used in the CKM-unitarity fits. Their explanation and discussion can be found, for example, 
in Ref. 25. 



6 Summary and Concluding Remarks 

We have investigated the impact of the current measurements of the time-dependent CP-asymmetry 
parameters S+~ and in the B^/B^ — > 7r + 7r~ decays, reported by the BELLE and BABAR collabo- 
rations, on the CKM parameters. The results of our analysis can be summarized as follows. 

In the first part of our analysis, we have compared the resulti ng wor ld average (J7J) of these measure- 
ments with the predictions of the specific dynamical approaches l ^ l ^* 4 ! in which the quantities \P C /T C \ 
and S c are estimated. We find that, within the SM, they do not provide a good fit of the data. In par- 
ticular, the estimates of \P C /T C \ and (5 C G23 based on the QCD factorization 1231 are off the mark by more 
than 3 sigma. This was shown in Fig. 0] for a large enough range of the angle a. The mismatch between 
the data and the QCD-factorization approach - - can also be studied by calculating the % 2 -function of 
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p 


0.10 - 


- 0.24 


n 


0.32 - 


- 0.40 




0.79 - 


- 0.86 


sin(2a) 


-0.44 - 


- +0.30 


sm(zp) 


0.69 - 


- 0.78 


sin(27) 


0.50 - 


- 0.96 


a 


(81- 


103)° 


P 


(21.9 -r 


25.5)° 


7 


(54 4 


75)° 


AM Bs 


(16.6 -7- 20.3) ps" 


\Pc/T c \ 


0.43 


5.3 


6 C 


(-112 - 


- -29)° 



Table 2. The 68% C.L. ranges for the CKM-Wolfenstein parameters, CP-violating phases, AMg g , \P C /T C \ and <5 C from the 
CKM-unitarity fits. 
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P 

Figure 11. Constraints in the p — fj plane from the five measurements as indicated. Note that the curve labelled as A.Mb s 
is obtained from its 95% C.L. lower limit 14.4 ps — . The fit contour corresponds to 95% C.L. and the dot shows the best-fit 
value. The black contour shows the impact of the GLSS lower bound cos(2(9) > 0.27 resulting from the C^T and B — * ittt 
branching ratios measurements and using the isospin symmetry. 



the unitarity triangle fit with the estimates by Buchalla and Safirl^for \P C /T C \ and S c as an input. The 
X 2 -minimum of the resulting fit is about 11 (compared to Xmin = 0-57, leaving these two parameters 
free), which corresponds to a probability of about 1%. Looking more closely to localize the source of 
this discrepancy, our fits show that it is the small value of the strong phase difference predicted in the 
QCD-factorization approach (—6° < S c < 24° at 68% C.L.) which should be compared to the fit of the 
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Figure 12. The 95% C.L. correlations between a — 7 and sin(2a) — sin(2/3) in the SM. The black contours show the impact 
of the GLSS lower bound cos(26>) > 0.27 resulting from the C^n and B — > tttt branching ratios measurements and using 
the isospin symmetry. 
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Figure 13. The x^; n - distribution as a function of the angle a from the unitarity fits of the CKM parameters in the SM 
including the current measurements of C^T" and S^^ ■ The dashed (solid) curve is obtained without (with) taking into 
account the lower bound cos(20) > 0.27. 



data obtained by leaving the two parameters free, yielding —64° < 6 C < —29°, which contributes mainly 
to the chi square and hence results in the poor quality of the fit. Since a similar inference also follows 
for the CP asymmetry Acp(K + tt ~) in the — > K + tt~ decay, where the current measurements yield ^ 
(-9.5 ± 2.9)%, compared to the QCD-factorization prediction (+5 ± 10)%, one must conclude, unless 
data change drastically, that this approach grossly underestim ates the strong interaction phases in the 
B — > 7T7T and B — > Kir decays. The competing pQCD approach E3I fares comparatively somewhat better 
but predicts a smaller value of |P C /T C | than is r equired by the current data. The central value of this 
quantity in the estimate by Keum and Sanda^HI \P C /T C \ = 0.23, is approximately a factor of two smaller 
than the lowest value of this quantity from the fit range 0.43 < \P C /T C \ < 1.35 at 68% C.L., with the 
best-fit value being \P C /T C \ = 0.77. The main message that comes out from this part of the analysis is that 
the QCD-penguin contributions are significantly stronger than most of their current phenomenological 
estimates, and it remains a theoretical challenge to understand this feature of the data. 

In the second, and larger part of this paper, we have addressed the question of how to interpret 
the measurements of S£~ and in terms of the CKM parameters in a model-independent way. We 
find that leaving the dynamical quantities \P C /T C \ and 5 C as free parameters, the CKM unitarity fits do 
not effectively constrain these parameters, yielding a very large range even at 68% C.L. As a result of 
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Figure 14. The Xm^-distributions as a function of [f c /T c | (left frame) and the strong phase difference <5 C (right frame) from 
the unitarity fits of the CKM parameters in the SM including the current measurements of C~%k an d S„„ . The dashed 
(solid) curve is obtained without (with) taking into account the lower bound cos(2f?) > 0.27. 
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Figure 15. Constraints in the p — fj plane from the assumed exact measurements of C~%^ = —0.46 and S^^ = —0.74 and 
the lower bound cos(20) > 0.27 (the shaded regions). The fit contour corresponds to the 95% C.L. unitarity fits without 
taking into account the C~^k and Stttt~ measurements. 



this, the measurements of and Cj~ have practically very little impact on the unitarity fits of the 
CKM-Wolfenstein parameters p and fj, and hence on the allowed values of the angles a and 7, unless 
these dynamical quantities are bounded. We have reviewed a number of proposals in the literature to 
put isospin-based bounds on the QCD-penguin contribution in the B — ► nn decays. Parameterizing it 
in terms of the angle 29, introduced by Grossman and Quinn, we find that the GLSS lower bound on 
cos(2#) is the strongest bound to date, which we have evaluated as cos(2f?) > —0.03 (propagating the 
errors on the input quantities) and as cos(2#) > 0.27 (for the central values). We have worked out the 
consequences of the lower bound cos(20) > 0.27 on the profile of the unitarity triangle and the angles a 
and 7. Including the isospin-based constraint, our best fit values yield: a = 92°, f3 = 24° and 7 = 64°, 
with the 68% C.L. ranges given in Table[5] The corresponding best-fit values of the dynamical parameters 
are found to be \P C /T C \ = 0.77 and S c = —43°, respectively, with their 68% C.L. ranges given in Eq. 155|) . 
With improved data expected in the near future, these ranges can be reduced significantly, leading to a 
precise determination of all three angles a, (3 and 7 of the unitarity triangle. 

Of course, at some stage, one has to take into account the isospin-breaking correction s. They o rigi- 
nate, in part, from the electroweak penguins which are estima ted to be numerically small I B | 45 | 46 | anc j 
this estimate can be put on model-independent gr 

oundsETEEl Moreover, they do not change the bounds 
obtained above for which the closure of the two triangles shown in Fig. [3 was used. However, the isospin- 
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breaking corrections may be significant from the t t° — r] — rf mixing I^Sa in the presence of which the two 
isospin triangles used in our analysis do not close leading to quadrilaterals. These latter corrections 
will have to be accounted for in the final determination of the angle a. As the current estimates of these 
corrections are model-dependent ESI we can not assign at present a quantitative weight to them. They 
will be better determined as and when the individual — ► 7r°7r° and i?° — + 7r°7r° branching ratios are 
measured to which we look forward in the future. 
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